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Abstract 

Using a stochastic generalization of the Hodgkin-Huxley model, we consider the 
influence of intrinsic channel noise on the synchronization between the spiking ac- 
tivity of the excitable membrane and an externally applied periodic signal. For 
small patches, i.e. when the channel noise dominates the excitable dynamics, we 
find the phenomenon of intrinsic coherence resonance. In this case, the relatively 
regular spiking behavior is practically independent of the applied external driving; 
therefore no synchronization occurs. Synchronization takes place, however, only for 
sufficiently large ion channel assemblies. The neuronal signal processing is thus likely 
rooted in the collective properties of optimally large assemblies of ion channels. 
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1 Introduction 



A fundamental question in neurophysiology concerns the limiting factors of 
the reliability of neuronal responses to given stimuli. In this article we fo- 
cus on a particular aspect of this complex issue: the impact of channel noise, 
which is generated by random gating dynamics of the ion channels in mem- 
brane patches of finite size. In particular, we investigate the effect of channel 
noise on the synchronization between the action potential produced by the 
cell membrane patch of finite size and the applied periodic stimulus. 
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The topic of synchronization, especially in biological systems, attracts ever 
growing interest, see Refs. [1,2,3,4]. Some prominent examples are the collec- 
tive flashing of fireflies [5], the synchronization between the respiratory and 
cardiac activity in human cardiorespiratory system [6] , and the signal process- 
ing in sensory systems [7,8]. Recently, the phenomena of frequency and phase 
synchronization have been explored in presence of ambient noise sources, both 
for a noisy (overdamped) relaxation dynamics [3,9,10] and for an oscillatory 
stochastic dynamics [3,11]. The theoretical research is further inspired by ex- 
perimental activities, and vice versa. In the context of theoretical modeling, 
the frequency and phase synchronization [12] has been observed, for example, 
in the integrate-and-fire model of the excitable dynamics driven by white noise 
and an externally applied stochastic spike train [13]. For an optimal dose of 
noise the mean firing rate of the driven neuron becomes locked by the mean 
frequency of the external spike train. 



2 The Hodgkin-Huxley model 



Our starting point is the well-established model of Hodgkin and Huxley [14]. 
The membrane patch of area S is considered as an electrical capacitor pos- 
sessing the specific area capacitance C. The membrane separates two ionic 
bath solutions (which in vivo correspond to the interior and the exterior of 
the excitable cell) with different concentrations of the ions of different sorts, 
mainly potassium, K + , sodium, Na + , and chloride, CP ions. The macroscopic 
concentration differences are kept approximately constant. In the cell, this 
task is accomplished by the ATP-driven ionic pumps. Furthermore, the ionic 
baths are on the average electrically neutral. However, due to the different 
ionic concentrations on the opposite sides of the semi-permeable membrane, 
the membrane becomes charged. As a consequence, an equilibrium transmem- 
brane electrical potential difference emerges. The lipid membrane creates an 
almost impenetrable barrier for the ions. However, they can flow across the 
membrane through special ion selective pores created by specialized membrane 
proteins - the ion channels [15]. The specific potassium, I K , and sodium, I Na , 
ion currents through the open ion channels are approximately proportional 
to the differences of the transmembrane potential V and the specific (for the 
particular sort of ions) equilibrium potentials, E K and E Na , respectively. The 
stochastically averaged, mean conductances, G^ a (m,h) and Gk(u), are, how- 
ever, strongly nonlinear functions of V . This nonlinearity emerges due to the 
gating dynamics (see below). There exists also the leakage current II, mainly 
due to the chloride ions. If the membrane is driven by the external current 
Iext{t), the sum of the specific ion currents and the capacitive current, Ic, must 
be equal to I ex t(t) as a consequence of the charge conservation. Therefore, the 
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equation for the transmembrane potential V(t) reads 



CjV + G K (n) (V - E K ) + G Na (m, h) (V - E Na ) 

+ G L (V - E L ) = 7 ext (0 . (1) 

For a squid giant axon, the parameters in Eq. (1) are Ej<s a = 50 mV, = 
— 77 mV, El = — 54.4mV, and C = 1 /iF/cm 2 . Furthermore, the leakage con- 
ductance is assumed to be constant, G L = 0.3mS/cm 2 . On the contrary, the 
sodium and potassium conductances are controlled by the voltage-dependent 
gating dynamics of single ion channels and are proportional to their respective 
numbers. These latter assumptions have been fully confirmed in the single- 
channel recordings by Neher, Sakmann and colleagues which indeed have 
proven that ion channels undergo the opening-closing stochastic gating dy- 
namics [16]. In the Ho dgkin- Huxley model, the opening of the potassium ion 
channel is governed by four identical activation gates characterized by the 
opening probability n. The channel is open when all four gates are open. In 
the case of sodium channel, the dynamics is governed by the three indepen- 
dent, identical fast activation gates (m) and an additional slow, so-termed 
inactivation gate (h). The independence of the gates implies that the proba- 
bility pK,Na of the occurrence of the conducting conformation is Pk — n 4 for a 
potassium channel and P^ a = m 3 h for a sodium channel, respectively. In the 
mean-field description, the macroscopic potassium and sodium conductances 
thus read: 

G K (n)=g^n\ G Na (m, h) = g^ x rr?h , (2) 

where g^ ax = 36mS/cm 2 and g^ x = 120mS/cm 2 denote the maximal con- 
ductances (when all channels are open). The two-state, open-closing dynamics 
of the gates is given by the voltage dependent opening and closing rates a x (V) 
and /3 X (V) (x = m, h, n), i.e. 

0.1Q/ + 40) 

amiV) = l-exp[-(V + 40)/10]' (3a) 

P m (V)=A exp[-(y + 65)/18], (3b) 

a h (V) = 0.07 exp[-(V + 65)/20], (3c) 

/3 h (V) = {1 + exp[-(V + Z^/IO]}- 1 , (3d) 

0.01 (V + 55) . , 

an{V) = l-exp[-(y + 55)/10] ' (3e) 

(3 n (V) = 0.125 exp[-(y + 65)/80] . (3f) 

Hence, the dynamics of the opening probabilities for the gates are given by: 

x = a x (V) (1 — x) — (3 X (V) x, x — m,h,n . (4) 
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The voltage equation (1), (2) and the rate equations of the gating dynam- 
ics (3), (4) define the original, purely deterministic Hodgkin-Huxley model 
[14] for the squid giant axon. Operating with the average number of open 
channels, disregards, however, the corresponding number fluctuations (the so- 
called channel noise [16,17]). Therefore, the Hodgkin-Huxley model is valid, 
strictly speaking, only within the limit of very large system size where these 
fluctuations can be neglected. The role of internal fluctuations in membrane 
patches of finite size cannot, however, be a priori neglected. In fact, the recent 
theoretical studies make it clear that the channel noise can be functionally 
important for the excitable dynamics [17,18,19,20]. 



3 Stochastic generalization of Hodgkin-Huxley model 



To account for the effect of channel noise, we use the model which presents 
a stochastic generalization of the Hodgkin-Huxley equations (1-4) due to Fox 
and Lu [21]. The dynamics of the gating variables is given by the following 
Langevin equation: 

x = a x (V) (1 - x) - /3 X (V) x + £ x (t), x = m,h,n, (5) 

with independent Gaussian white noise sources £ x (t) of vanishing mean which 
take into account the fluctuations of the number of open gates. The noise 
strengths depend on the membrane voltage. The noise correlations assume 
the following form for an excitable membrane patch with N^a sodium and 
potassium ion channels [21]: 

VVNa ("m + Pm) 
^VNa \Ot h + Ph) 

To ensure the confinement of the gating variables between (all gates are 
closed) and 1 (all gates are open) we implemented numerically the reflecting 
boundaries at and 1. With the assumption of homogeneous ion channel 
densities, p^ a = 6 • 10 13 m~ 2 = 60 /im~ 2 and pk = 1-8 • 10 13 m~ 2 = 18 /iin~ 2 , 
the ion channel numbers are given by: 

iV Na = PNaS, N K = p K S . (7) 

With decreasing patch size S we observe from our numerical simulations of 
the Langevin equations (5-6) the more and more increasing spiking activity 
due to the internal channel noise, see Fig. I. 
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Fig. 1. Three simulated realizations of voltage spike trains for different patch sizes: 
(a) 1 /im 2 , (b) 16 /im 2 , and (c) 128 /im 2 . With decreasing patch size (bottom-to-top) 
the fluctuations of the number of open ion channels dominate the dynamics; there- 
fore more action potentials are produced, whereas for very large patch sizes the 
intrinsic channel noise strength is not capable to produce corresponding voltage 
spikes. 

4 Synchronization in excitable membranes 



4-1 Channel noise induced spiking activity: Intrinsic Coherence Resonance 



We investigate numerically the influence of channel noise in the absence of ex- 
ternal stimulation. For small internal noise strengths, i.e. large patch sizes, 
the occurrence of action potentials is very rare in the autonomous, non- 
driven regime: J cxt = 0, see Fig. 1. The distribution of interspike intervals 
{Tj := ij+i — ti} is widely spread with longer intervals becoming less and less 
probable, see the solid line in Fig. 2(a). With increasing channel noise the fir- 
ing rate increases as well, and the interspike interval time-histogram depicts 
a more distinct peak structure, Fig. 2(b). For small membrane patches, the 
fluctuations of the number of open ion channels dominate the excitable dy- 
namics and after a firing event the membrane patch tends to fire immediately 
again, but after some refractory time which is necessary to achieve the resting 
potential after hyperpolarization. The above elucidation suggests that with 
decreasing patch size the spiking activity becomes more regular. A quantity 
that measures this regularity is the so-called coefficient of variation (CV), or 
the relative dispersion of the interspike interval distribution, i.e. 

which involves the mean interspike interval (T) := lini/v^oo — U)/N 

and the mean squared interval (T 2 ) := lirn ; v- >00 J2(U+i — U) 2 /N. For a fully 
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Fig. 2. The interspike interval histograms (ISIH) for the undriven case, I ex t = 0, 
(solid line) and the driven case, I ex t = sin(0.3 • t), (dotted line) for two patch sizes: 
(a) 16 /on 2 , and (b) 1 /an 2 . The driving independent single-peak-structure becomes 
more and more accentuated for decreasing patch size, i.e. with increasing strength 
of the channel noise. The refractory period T r acts as lower bound to the inter- 
spike intervals. It varies with the patch area and becomes smaller for smaller patch 
sizes, cf. T r fa 15 ms in part (a) vs. T r < 10 ms in part (b). The driving-induced 
peaks appear for sufficiently small internal noise at multiples of the driving period 
^driving ~ 21 ms, cf. part (a). In case of dominating channel noise, the influence of 
driving on the ISIH is barely distinguishable, cf. part (b). 

uncorrelated sequence of spikes, which corresponds to the Poissonian spike 
train, the coefficient of variation would assume the value CV = 1. For a more 
ordered spike train, the coefficient of variation assumes values less then one, 
CV < 1. For a purely deterministic signal, the CV equals to zero. 

Fig. 3 depicts the coefficient of variation (8) and the firing rate k = 1/(T). 
While the firing rate decreases monotonically with respect to the patch area, 
the CV reveals the novel phenomenon of intrinsic coherence resonance [19]. 
At an optimal dose of internal noise, i.e. an optimal size of the cell membrane 
patch near S = 1 /mi 2 , the CV exhibits a minimum, where the spike train 
becomes distinctly more ordered. With decreasing noise level (large S) this 
feature becomes increasingly diminished, and the CV increases. 



4-2 Sub-threshold sinusoidal driving 



Next, we switch on an external, sub-threshold sinusoidal driving, i.e. 

7 ext := Asm(nt). (9) 

with A < A th , where A th is the threshold amplitude of driving in the de- 
terministic limit, S — > oo. For Q = 0.3 ms -1 in Fig. 2, A t h ~ 1.6 /zA/cm 2 . 
Interestingly enough the distribution of interspike intervals is not affected 
much for small patch sizes near S ~ 1 /im 2 , cf. Fig. 2(b), where the solid and 
dotted curves are barely distinguishable; see also Fig. 3(a) at 5 « 1 /im 2 . In 
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Fig. 3. (a) The dependence of the coefficient of variation CV in (8) versus patch 
sizes, (b) Firing rates k = 1/{T) for the undriven (solid line) and driven case (dotted 
line, I ex t = sin(0.3 • t)) are depicted versus the patch areas. While the firing rate 
is monotonically decreasing with increasing patch size the coefficient of variation 
evidences the phenomenon of intrinsic coherence resonance. 

this case, the spiking activity possesses an internal rhythm which withstands 
the external disturbances. For larger patch sizes the internal noise decreases 
and the periodic driving induces peaks at multiples of the driving period, see 
Fig. 2(a). In this regime the external driving imposes a more ordered spiking 
activity which is characterized by a reduced CV, see Fig. 3(a). 



4-3 The concept of stochastic phase processes: frequency synchronization 



In order to explore the stochastic phase and and stochastic frequency synchro- 
nization between the membrane spiking activity and the driving current, one 
has to invoke the concept of a phase process linked to the stochastic voltage 
signal V(t). The Hilbert phase H (t) constitutes one such appropriate phase 
definition. The corresponding approach is based on the Hilbert transform, it 
has been originally introduced by Gabor [22,23]. The signal is extended into 
the complex plane with the imaginary axis given by the Hilbert transform of 
signal, V R , cf. Fig. 4: 

z v (t) :=V(t)+iV u (t), (10) 

with 

V«(t) = H[V](t) = ±P J ^dr. (11) 

71 J t — T 

— 00 

The integral in (11) is evaluated in the sense of the Cauchy principal value 
(P). The Hilbert transform is then used to define the corresponding stochastic 
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Fig. 4. Stochastic (Hilbert) phase process: (a) realization of the voltage signal (solid 
line) and the corresponding Hilbert transform (dotted line); (b) The trajectory 
in {V,V H }-space is plotted for the time interval in (a) which is determined by 
two dashed vertical lines at t = 40 ms and t = 80 ms, respectively, containing a 
single firing event. The {V(t), V H (i)}-space trajectory undergoes a full 27r-revolution 
during such a single firing event. 

process for the Hilbert phase H (t), i.e. 

v H {t) 



i H (t) = arctan 



V(t) 



(12) 



Here, the phase should be understood as continuously growing function of 
time. In Fig. 4(a) a realization of the voltage signal and the related Hilbert 
transform are plotted. Each repeated firing event corresponds to a cycle in the 
{V, V H }- space, see 4(b)); it thus adds an 27r-increment to the phase of the 
voltage signal. The corresponding Hilbert phase frequency wh is then given by 



lim 

t—>oo 



<t> K (t) 



(13) 



An alternative phase definition uses the feature of the point process that is 
generated by the spikes of the voltage signal V(t). These spikes define marker 
events at times {ti} with each of them pinpointing the completion of a cy- 
cle. The phase increase between two subsequent marker events is then given 
exactly by 2n. Moreover, by a linear interpolation it is possible to define the 
instantaneous phase R (t); i.e., 



R (t) = 2tt 



t-U 

ti+l ~ ti 



+ 2vri , (U < t < t i+1 ) , 



(14) 



where the times ti are fixed by the marker events. The average frequency of 
the process (f) K (t) is then given by: 



Wr = lim 



R W 



t^oo t 



2it lim 



N(t) 



t^oo t 



(15) 
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where N(t) denotes the total number of spikes within the time interval [0,t]. 
The index R alludes to the so-called Rice frequency [3,11,24]. 

Since each action potential adds, on the one hand, the 27r-increment to the 
Hilbert-phase and, on the other hand, defines the corresponding single marker 
event, the Hilbert frequency equals the Rice frequency, although the two phase 
definitions are different. 



4-4 Imperfect synchronization for a periodically driven membrane patch 



0.4 
0.3 



3 0.2 



0.1 
0.0 



vja) 




^ X 











1 10 10 2 10 3 10 4 




S [fim 



tt/2 71 3tt/2 2tt 
ttt [mod 2tt] 



Fig. 5. (a) The dependence of Rice frequency (15) on the patch size is plotted for 
three sub-threshold amplitudes A of the sinusoidal stimulus I ex t(t) = ^4sin(0.2 • t): 
A = 0.0 uA/cm 2 (short-dashed line), A = 2.02 uA/cm 2 (dotted line), A = 2.05 
/xA/cm 2 (long-dashed line), and for a super-threshold amplitude A = 2.2 fiA/cm 2 
(dotted-dashed line). The threshold amplitude is A t h « 2.1 uAjcm 2 . The level line 
displays the angular driving frequency Q = 0.2 ms _1 . For subthreshold signals the 
firing rate vanishes for large areas and the Rice frequency approaches zero, while 
for a super-threshold amplitude the driving frequency value is attained, (b) The 
phase probability density p{4>dr) of spiking events versus the phase <pdr = &t of 
the driving signal is plotted for A = 2.05 fiA/cm 2 and four different patch areas: 
S = 4 am 2 (solid line); S = 16 fim 2 (short-dashed line); S = 64 ^m 2 (long-dashed 
line); S = 256 fim 2 (dotted-dashed line). 

The dependence of the mean frequency wr on the patch size is depicted in Fig. 
5(a) for driving signals with frequency Q = 0.2 ms _1 and for different driving 
amplitudes A. In the strong noise regime, the spiking activity is determined by 
the internal noise and is independent of the driving frequency and amplitude. 
For intermediate-to-small noise strengths the spike occurrences are locked to 
multiples of the driving period, cf. Fig. 2 (a). 



This synchronization behavior is also reestablished in the phase probability 
density p{<pdr) of spiking events in Fig.5(b). The latter quantity is defined 
with respect to the phase (pdr = ^ of the driving signal within the driving 
period, i.e. < (fidr < 27r. For small patch sizes, i.e. for a large channel noise, 
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the probability density p{<pdr) is relatively flat. When the channel noise di- 
minishes a peak structure becomes more and more pronounced. It worth to 
notice that spikes occur most frequently before the driving current reaches its 
maximum at fit = n/2. The occurrence of the spike thus forecasts that the 
signal's maximum will be reached soon. This effect is similar to the effect of 
anticipated synchronization between master and slave neurons [25] . Moreover, 
a very interesting new effect is noticeable in Fig. 5(b): the phase lag between 
the maximum of p((fidr) and the signal maximum at <pdr = vr/2 increases upon 
increasing the channel noise strength. In other words, the smaller membrane 
patch fires most probably more in advance than the larger patch, although the 
probability distribution of the spikes occurrence is increasingly flattened. This 
prominent effect can be explained as follows. The increase of noise strength 
certainly helps to overcome the firing threshold. Therefore, the system tends to 
fire more and more in advance, before the maximum of signal is reached. How- 
ever, the noise increases also the uncertainty of firing events and, therefore, 
the probability distribution becomes flat. 



Even though the observed locking behavior presents clearly some sort of syn- 
chronization, no perfect frequency synchronization - characterized by a per- 
fect locking of the mean frequency (13) to the external driving frequency - 
can be detected, see Fig. 5(a). While for large, but still sub-threshold stim- 
uli A < A th pa 2.1 //A/cm 2 a plateau is formed over a wide range of noise 
strengths, the mean frequency, which corresponds to this plateau, doesn't 
match the driving frequency. This frequency mismatch happens due to the 
multimodal structure of ISIH, which is caused by the locking of the firing oc- 
currences to the external force period in ratios different from 1 : 1. A similar 
phenomenon of imperfect synchronization has also been found in the human 
cardiorespiratory activity [6]. The plateau corresponds to the situation where 
nearly every maximum of the driving force produces a spike. Nevertheless, 
some driving maxima produce no spikes and therefore the mean frequency is 
smaller than the driving frequency. We point out, however, that for a signal 
represented by a stochastic spike train, like in Ref. [13], a perfect synchroniza- 
tion should be observed. The reason is that a random driving with the mean 
frequency (fl) will not produce additional maxima in the ISIH at the multi- 
plies of the mean driving frequency. The occurrence of such maxima thereby 
destroys the perfect synchronization in the periodic case. 



The situation changes with a super-threshold stimuli, see dashed-dotted line 
in Fig. 5(a): here the increasing strength of the channel noise with decreasing 
size S destroys the perfect frequency synchronization. This is caused by the 
skipping of some firing events, or by the generation of additional spikes. 
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5 Conclusions 



In conclusion, we have investigated the synchronization in a noisy general- 
ization of the Hodgkin-Huxley model, which incorporates the spontaneous 
fluctuations of the membrane conductivity due to the individual ion channel 
dynamics - the so-called channel noise. In the absence of an externally ap- 
plied stimulus the excitable membrane patch exhibits a noise-induced rhyth- 
mic spiking activity at an optimal patch area. This effect can be regarded as 
the intrinsic coherence resonance [19,20]. In the presence of external periodic 
driving, the interspike interval histograms clearly demonstrate a synchroniza- 
tion of firing events with the external driving for large and intermediate patch 
sizes, while for small membrane patches the channel noise reigns the spik- 
ing dynamics. Due to the excitation of higher harmonics, a perfect frequency 
synchronization could not be observed. Instead, the Rice frequency, which is 
defined by the occurrence frequency of action potentials, exhibits a plateau-like 
structure for subthreshold signals: over a rather wide range of channel noise 
strength the Rice frequency can be fixed at a value smaller than the driving 
frequency. Furthermore, it was shown that the Hilbert frequency of the voltage 
signal equals to the Rice frequency of the spike process and, therefore, this 
average Rice frequency also reveals an imperfect synchronization behavior. 
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ter on the nanoscale, Project A10. 
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